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SmWARY 


The shook-expansion method of NACA TN No, llli3 was used 
to detemine the crinoioal aerodynamic characteristics of 
two-dimensional supersonic airfoils. A discussion is given 
of the effect of thickness ratio, free-stream (Jach number, 
angle of attack, camber, thickness distribution, and aileron 
deflection. The calculations indicated that the minimum 
drag of supersonic airfoils is obtained when the maximum 
thickness is behind the O. 5 O chord. The center of pressure 
obtained for a symmetrical supersonic airfoil was found to 
be ahead of the O .50 chord. 


INTRCDDCTICH 


The characteristics of thin airfoils moving at supersonic 
speeds are determined in reference 1 by Ackeret's thin-airfoil 
theory. In this method, the local static pressure is assumed 
to vary linearly with the angle between the free-stream 
direction and the local airfoil surface. This assumption 
precludes any effect of camber on lift and locates the center 
of pressure of an uncambered airfoil at the nidchord. 

The relations for flow across shock waves are presented 
in reference 2. Reference 3 ocanbines these shock equations 
with Meyer’s expansion equations (see reference I4,) and presents 
a graphical way of calculating a second-order approximation 
to the characteristics of thin sborp-nose airfoils at super- 
sonic speeds. The present paper uses the shock-expansion 
method of reference 5 "to determine sesme interesting effects 
of thAckness ratio, free-stream Mach number, angXe of attack, 
camber, thickness distribution, and aileron deflection, 
Swept-back wings are not considered herein as they are of 
sufficient interest to justify a separate report. 
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SYMBOLS 

M Manh niamber 

p static pressure 

q dynamic nreesure 

a anp;le of attack of airfoil 

c section lift coefficient 

V 

c^ section drag coefficient 

coefficient of section pitching moment about airfoil 
•5o raidohord 

t airfoil thickness 

0 airfoil chord 

o„ aileron chord 

a 

6 aileron deflection 

Subscripts: 

0 free stream 

1 unner leading surface 

2 upper trailing surface 

3 lower leading surface 

if. lower trailing surface 

mar maximum 



DISCUSSION 


Airfoils experience an increased, oressurf: draft at high 
Mach numbers. This drag increase can be minimized by the use 
of thin airfoils with sharp leading edge^s. For this reason, 
and also because theoretical calculations are more accurate 
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for this type of airfoil, the wing sections considered herein 
are limited to fairly thin, sharp-nose airfoils having little 
surface curvature. 

Center of nress’Jtre.- Reference 5 has given the pressure 
distribution at = ^4. the diamond -shape (syrametrio double- 
wedge) airfoil shown in figure 1 . The pressiure coefficients 
around this airfoil are given in, table I, as follows: 

TABLE I.- IRESS^E COEFFICIENTS 


Location of surface 
(see fig. 1) 

Pressure coefficient, 
P-Po 

9 o 

Upper leading 

-0.0169 

Lower leading 

.024.16 

Upper trailing 

-.0308 

Lower trai ling 

.0186 

( 


A study of this table brings out many interesting points 
on the characteristics of supersonic airfoils. The lift of 
the airfoil is proportional to the difference in the pressure 
coefficients on the upper and lower surfaces of the airfoil. 

For the leading half of the airfoil of figure 1 , this difference 
is 


O.Ci4-l6 + 0.0169 = 0.0585 

For the trailing half of the airfoil the difference is 

0.0308 + 0.0188 = O.CI4.96 

The significant result to be noted is that the front of the 
airfoil is carrying more than half the lift and the center 
of pressure is found to be at the i^6-percent-chord position. 

As the supersonic Mach number is decreased toward 1 . 0 , the 
center of pressure of thin diamond-shape airfoils approaches 
the 50“P3rcent-chord position. As the supersonic Mach number 
is increased, the center of pressure of thicker airfoils moves 
forward. Figure 2 shows the variation of the center of pressure 
of a 5“percent-thick diamond-shape airfoil with angle of attack 
and Mach number. The actual shift in center of pressure depends 
on the airfoil shape. For example, a wedge airfoil having its 




k 


me A TN No. 1179 


raaxinsum thickness at the trailing edge v;ill have its center of 
pressure at midchord at all supersonic speeds provided that the 
angle of attack is not sufficient to cause the shock wave to 
separate from the leading edge. 

Drag coefficient. - In order to determine the actual 
variation of ura^, soine allowanoe must be made for the viscous 
skin friction. Tlteoreti cal work (reference 2) suggests that the 
viscous drag coefficient decreases with increasing Re^molds 
number at supersonic speeds. However, reference 5 shows that the 
viscous dr'^g coefficient Is almost unaffected by changes in Mach 
number or Reynolds ntamber if the Roynoldb number is, very high. 

For this reason the skin-friction drag oo"ffioient based on wing 
area is considered herein to be constant: at 0.0060, which is 
approximately the value obtained for reefent highly polished 
jet-propelled fighters; however, the results are plotted in 
such a way that the reader oan easily adjust the curves to 
corrr sDcnd wi.th the conditions in which he is Interested. If 
new information indicates that the variation of friction with 
?J?ach numher is appreciable, the curves may be raised or lowered 
by the amount of the variation. 

Fi?\rre 3 shows the- variation of the drag ooeffioient of 
diamond-shape airfoils with thickness'^ratlo and free-stream 
Tiaoh nuDiber. The eranh shows that the drag coefficients are 
very high near a l?aoh number of 1.0 and the main part of 
the drag then is pressure drag. At high iiach numbers the shock 
drag has decreased in importance relative to skin friction. 

It is seen that the Aokeret method (reference 1) gives almost 
the same trends us the shook-expension method except that it 
shows a less rapid variation of drag with thiokness at high 
T'iaoh numbers than the present shock-expsuaslon method. Ackeret's 
method predlots lower pressure increases on the leading portion 
of the airfoil and higl.er pressure decreases on the trailing 
portion than the method of combining the ’shock and expansion 
relations. 

Ihe effect of the location of the ;\axlmum thickness on 
the drag coefficient for a 5-psrcent-thick airfoil is shewn 
in figiire ij.. This figure indicates that the optimum location 
of maximum thiokness is close to the midchord position for 
fairly low supersonic speeds. At Mach n’Jmbers of 6 and above, 
however, the optimum location seems to ba near the trailing 
edge. This condition is very different from that predicted by 
Ackeret's method whore, by the nature of .his assumptions, the 
pressure distribution is symmetrioal. In connection with airfoils 
designed for an angle of at’back other than sero it must be kept 
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in mind that the optimm location of thickness on the upoer and 
lower surfaces will he different. 

Figure 5 shews the effect of location of maximum thickness on 
the drag of IQ-peroent- thick airfoils. These airfoils show an even 
greater variation from Aokeret’s prediction than the 5 ~?ei*< 3 ent- 
thiok airfoils. This figure indicates tliat the optimum position 
for maxim’im thickness at a Mach number of B is at the trailing * 
edge. Actually, the flow at the trailing edge will break down 
for this condition and form a turbulent wake of somevd'.at higher 
pressure than that calculated. This breakdown of the theory for 
the airfoils \-^th hliint trailing edges should make the experimental 
drag less than the theoretical drag for this extreme condition. 

Lift coeff iolent . - In figure 6 the slope of the lift curve is 
clotted as a ilinotion of frec-stream Mach number for three different 
airfoil thicknesses. At a Mach, number slightly above 1 th? slope* of 
tho lift curve is apuroximately equal to the Gubaonio slope. The 
slope of the ■ lift curve drops rapidly, however, with increasing Mach 
number. The trend of the curve is similar to the drag curve of 
figure 3 except that the drag curve was displaced upward a constant 
am.ount by the skin friction. The thickness ratio of the airfoil 
seems to have little effect on the slope of the lift curve except 
at high supersonic Mach numbers (above Ij.) . The thick sections 
then have the highest slope. As might be expected, the calculation 
obtained by the shock-expansion uiethod for airfoils of zero' thick- 
ness give results identical wnth those given by Ackeret’s method. 

Li.ft-drag ratio. - The effect of angle of attack and Mach 
number on the characteristics of a thick double-wedge 

airfoil is shown in figure 7« This figure shows that the lift 
coefficient increases almost linearly with angle of attack, and 
the drag, coefficient increases with angle of attack in a marner 
similar to tlie variation of total drag coefficient for a complete 
subsonio wing. 

The maximum lift-drag ratio oj/Cji a-t a Mach number of 2 is 
approximately 6.92 at c, wO.lo and a = I 4 P. At a Mach number of I 4 . 
the maximim ratio is 5«52 at Oj = 0«09i4. and a = 5°* The maximum 
lift-drag ratio decreases with increasing Mach number, and the lift 
coefficient for maximum lift-drag ratio also decreases with in- 
creasing l!aoh number. These trends are somewhat different from 
those of reference 1 in which skin friction was neglected. In that 
case the maximum lift-drag ratio was independent of Mach number. 
Bussmann in reference 6 calculates the Mach number at which the 
lift-drag ratio is maximum for certain airfoil thicknessos. 

Camber . - Table 1 has shown that the leading portion of an 
airfoil -bends to carry more than its share of the lift and, hence, 
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it can be expected that camber, ishlch decreases the angle of attack 
of this impcrtaut part of the airfoil, will decraase •fehc lift. In 
figure 8, calculations show that the addition of 1.25 percent comber 
changes the characteristics piTren in figure 7 for i'ach number ij.. 

The addition of coLiber increases the dfas coefficient at zero angle 
of attack but decreases the drag ooeffioient at. high angle of attack. 
The angle of zero lift is shifted slightly but the slope of the lift 
curve remains essentially the same. It is significant that the 
shift obtained as a result of camber at suncrsonic speeds is opposite 
in direction to the shift obtained at subsonic speeds. This result 
indicates thuit a cambered wing will experience a change in the angle 
for zero lift in ao celt rating through the speed of sound. The 
addition of camber decreases the. lift-drag ratio at low angles of 
attack o.rid increases it at high angles -of attack. 

A comparison of the angle of uttaok for zero lift as calculated 
by the AolTeret method (reference 1), as calculated by the shock- 
expansion method of reference determined experimentally 

by referonce 7 is pro sonted in table II. Thus, 

TABLE II.- ANGLE OF /ATTACK FOR EEHO LIFT AT L'ACF K’JJiBSR 2.15 


J 

Airfoil shape 

Angle of zero lift, deg | 

_ ■ _J 

Ackeret 

theory 

'reference 1) 

■ - r "■ f 

Shook- i 

expansion 1 

method ! Experimental 

(reference 5) [ (reference 7) 

G.U.5 

0 

1,68 

^ 1 
1.89 


0 


.45 

6^10' 

kPio^ 

■ 

0 

.28 

- - 

.26 
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The shock-expansion method (reference 5) shows better agree- 
ment with experimerit tlr^an bhe Ackeret method as far as the effect 

of camber is concorned. Theoretically, a slight error might be 

expected for bhe G.U.3 airfoil because of tVie inter feronco of 
shock and expansion waves due to the continuous surface curvature; 
howp-n-er, the error seems to be of very low magnitude. 

In connection with table II it might be well to 
point out that higher order terms have been developed 
by other authors which, when added to the Ackeret 
theory, predict soma effect of camber on lift, 

Ai leron s. - In order to dsmor. etrute tho effect of ail-rons on 
super scnic airfoils, results are giv^n in fig>ire 9 to show the 
aileron effectiveness factor for' a 5~P®rcent-thiok diamond- shape 
airfoil at a Vaoh nmber of i; us a function of the ratio of aileron 
chord to airfoil chord. For comparison, a curve for the seme 
airfoil is given for the subsonic ( incccipressible-flow) condition 
(Y — 0) , us well as one obtained by Aokeret's method. The 

do7/d6 

aileron effectiveness factors — ie are shown to be much 

lower at suoersonio speeds than at speeds where there ai*e no 
compressibility effects. These v,alues are also slightly lower than 
those obtained by the Ackerct method. For instance, a 20-percent- 
chord aileron at Yach rtimber has anproximately tho same effective- 
mess factor as a l-percont-chord aileron at low subsonic speeds. 

This result indicates that the helix tingle of on airplane in 
rolling will probably be muck less at supersonic speeds than the 
helix angle described by airplanes at low subsonic speeds. At 
these speeds bhe aileron influences the wing ahead of it, but at 
sunersonic speeds the aileron alone is affected. In. fact, the 
wir-.g at supersonic speeds loaves the aileron in low-density air, and 
thus the cffectiviness is decreased to values lower than might 
otherwise be expected. The actual effoctivensss varies for 
different" airfoil shapes and thickness ratios. Fig’ore 10 gives 
a c\irve of the variation of section lift ooefficienb due to aileron 
deflection and figure 11 gives the slope of this curve; a slight 
increase in aileron ci'fectivenecs with increasing deflection is 
indicated. 

The section pitching-moment coefficient e.bout the midchord 
Cjji as a function of aileron dof lection 6 is . shown in figure 12. 

In figure 15 the slope of this curve dCjjj /d6 is plotted against 

.5q' 

aileron deflection and shoi^s an increased rate of change of pitching- 
moment coefficient with increasing aileron deflection-, The 
pitching-moment coefficient about tho midchord is plotted in figure lit 
as a f’XQction of angle of attack for constant aileron deflection. 

The absolute valxje of the pitching-moment coefficient is shown to 
decrease with increasing angle of attack for a downward aileron 
deflection of 10*^ at F Ii. 
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CONCLUSIONS 


Calculations made by the shock-exoansion method to determine 
the aerodynamic characteristics of supersonic airfoils indicated 
the following conclusions ; 

1* Unoombered double-wedge airfoils having the maximum thick- 
ness at midchord -will have thoir centers of nressure ahead of 
the- 50 -percGnb chord at high supersonic speeds, but this center 
of pressure will approach midohord as the Mach ntmiber is lowered 
toward 1.0. Airfoils having their maximum thiokness near the 
trailing edge will have the center of pressure near the midohord 
at all supersonic speeds provided the angle of attack is not 
suffioient to oause the shock wave to separate from the leading edge. 

2. Tho pressure-drag coefficient and the lift coefficient for 
the same angle of attack decrease in a similar manner with increasing 
?*-aoh number, and thus their ratio is essentially constant with Mach 
number. The addition of a constant skin-friction drag coefficient 
results in a decrease in the lift-drag ratio with increasing 
Mach number, 

5* The optimum location of maximutn thickness for a given 
+hiokncss ratio to give minimum drag depends on airfoil shape and 
frec.-sfcream f'ach munb<.-;r. For double-wedge airfoils the optimum 
position of maximum thickness is near the trailing edge at very 
high Mach numbers 3 however, the ootimiim position approaches the 
miiiobord as bhs speed is decreased toward a Maoh number of 1.0. 

ij.. The aileron effeotiveness factor is lower when estimated by 
the shock-expansion method than when estimated by the Ackerct method. 

5* A comparison of theoretical and experimental values of the 
angle of zero lift suggestB" that the present method of calculation 
is a oonsiderably closer approximation than Ackerct' s method. 


Langley Memorial Aeronautical Laboratory- 

National Advisory Committee for Aeronautics 
Langley Field, Ya., July 18, 19i;6 
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Flgnr* 2»- Bosition of eontor of Tptm»tmv for an uBoaiA>erod 

airfoil* 


5- pareant-thlek 
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Fig. 3 
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Fig. 5 



Looatloa of MTium ‘UiioknasSi pereon't chord 

Figure 5*«> Effeot of looation of aaxiBUin thiokneaa on drag of 
10 “ p«roent>thlok airfoil. 



Slope of the lift enrfe per dag, de,/da 


Fig. 6 
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rree*stre«B Mteh auaberf 

Flgnre 6«- Slope of lift eu rre of oneuibeared thin airfoil at enpereonle epeede* 

t^Bx aidebordi s 0* 
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Fig. 7 


0 ^ 

c- ~r Airfoil contour 

Oj/«d 



Anglo of ottook, oC. , deg 


Figure 7«- Effect of angle of attack and Uaoh nuaber on charaeteristioa of 
5 - peroant>thlok unoaabered airfoil, t— »• et eidohoM} skin friction 

drag coefficient^ 0.0060. 


’ag ratiOj 



Section dx«g coeflleieBt 


Fig. 8 
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i — Airfoil with 0 eoaber 

(from figure 7) 

Airfoil with 1,25 percent 

CMiber 



Angle of attack, oC , deg 


Figure S.- Effect of angle of attack and canber on characterietles 
of 5 - percant-tblok airfoil at 1^ : A, 
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Subsoale 

(M^ « 0) 


Shoek>«zpansioa 


^ ^ . NATIONAL ADVISORY 

Aileron ohordAing chord, e./o COMMITTEE FOR AERONAUTICS 

Figure 9«* Aileron effeotlreness as a function of the ratio of 
aileron chord to airfoil chord for 5 * peroent>thiok unoaolbered 
airfoil at Ho s 4 , t ^^ at aidohord, 5 percent. 
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Angla of aileron deflection, cT, deg 

Figure 12*« Effect of aileron deflection on pitohing-moaent 
coefficient of airfoil section, tj^^x nidohord, 

5 percent; camber, Of cc m cPf • 4 * 


Rate of change of airfoil pltohlng>noaent 
coefficient with aileron deflection, do^/d^ 
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Fig. 13 



Angle of aileron deflection, , deg 

Figure 13 •- Effect of aileron deflection on elope of airfoil pitching 
Donent-ooeffloient curve for 20 percent chord aileron at Mq 3 4 
tgjg^ at nidohord, 5 percent; caaher, 0; cc s cP, 





